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NATIONAL ADVISORY "COMMITTEE FOR AERONAUT ICS 



TECHNICAL MEMORANDUM NO. 1018 



ON TEE ..SYMMETRICAL POTENTIAL FLOW OF COMPRESSIBLE 
FLUID PAST A CIRCULAR CYLINDER IN THE TUNNEL 
IN THE SUBCRITICAL ZONE* 
. By Ernst. Lamia 



The two-dimensional symmetrical potential flow of 
compressible fluid past a circular cylinder placed in the 
center line of a straight tunnel is analyzed in second 
approximation according to the Jansen-Rayl eigh method. 
The departure of the profile from the exact circular 
shape can be kept to the same magnitude as for the incom- 
pressible flow. The velocities in the narrowest section 
of the tunnel wall and at the profile edge are discussed 
in detail. 

In the treatment of the two-dimensional potential 
flow of an incompressible fluid past a circular cylinder 
placed in a tunnel with straight walls, the . ref 1 ect ion 
method is resorted to; the doublet which may be envisaged 
in the center is reflected at the tunnel walls and the 
effects of all the image points added. The then obtained 
solution is exact; the profile washed, tj the air stream 
is, to be sure, no longer an exact circle, but the differ- 
ence is ordinarily quite small. This solution for the 
incompressible flow is designated as the first approxima- 
tion; the flow of fluid past a circular cylinder in the 
tunnel with allowance for the compressibility effect, 
discussed in the present article, forms the second approx- 
imation. The study is based on a method originally pro- 
posed by 0. Janzen (reference 1) and soon afterward quoted 
by Lord Rayleigh (reference 2)**. 



*"Die symmetr ische Potent ialstromung eines Kompressibeln 
Gases urn einen Er ei szylinder im Kanal im unt erkr it ischen 
Gebiet." Luf tf ahrt for schung , vol. 17, no. 10, October 28, 
1940, pp. 329-331. 

♦♦Potential cp and stream function \j/ for the circular 
cylinder (and the sphere) in the free gas space have been 
computed up to the fourth approximation by this method 
for the first time by the present writer. 
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Ihe considerations prbee*d from the stream function 
\J/ with p denoting density, p pressure w speed, a 
corresponding velocity of sound, and ao the velocity of 
sound in the medium at rest. In absence of mass forces 
Euler's .equations for stationary, parallel flow afford 



2 " J P 



const 



(1) 



With the introduction of the stream function i|/ the 
equation of continuity div(pw) =0 is satisfied by the 
formula 



p 9y 



p S x 



(3) 



where U is the flow velocity p^ the density at infin- 
ity (undisturbed flow). The center line of the tunnel is 
the X-axis. By reason of the nonvorticity of the entire 
motion rot z w = 0 which, after posting (2) , gives 



A\J/ = grad , \J/ grad log p 
since ordinarily • dp = a 8 dp it follows from (l) that 

|- grad •if + a s grad log p = 0 



(3) 



whence, according to (3), 



grad w 3 grad \]/ - 0 



(4) 



If the regarded processes are adiabatic we get 



with K = ~ 

C T 



K -1 3 

W s 

2 ~ 



(5) 
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ttext we assume with Janzen and Rayleigh that quantity 
ii can he expanded in. a series, which progresses in powers 
of ratio a from flow velocity U to sonic velocity a Q . 
(In this instance q4/= u / a o *s substantially in agreement 
with the Mach number for undisturbed flow li = u/a^. ) 
The structure of equations (4) and (5) is such that only 
even powers of a occur. Hence we put 

a » U/a Q 



(6) 



n=0 

CO 

W : 2/U 2 = E = E^a 2 * 
n=0 

where ^ n and K n are" no longer dependent upon U and a ( 

Replacement of a s in (4) by (5) followed by posting 
(S) , the expansion in powers of a on both sides in com- 
parison of the coefficients of equal powers on both sides 
affords _ 



A\|> 0 = 0 



n 

AMh+1 = | Y (-gra# m grad E n .. m + (ft-l)E n _ m A^ m ) 
m=0 



(7) 



from which the a|/ can be deduced consecutively. The pre- 
diction of K n from \{/ n by means af (2) and (5) and the 
adiabatic equation f o Hows Yrf"^" 1 ^ 




Because of the appearance of density p the solution 
of the velocities from is a little more complicated 

than that from cp. But the use of ^ merits, because of 
the direct relationship with the profile form, preference 
in many instances. For the second approximation, that is, 
the prediction of ' ty x only quantity So is employed, 
for which (3), (6), and (7) afford 
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o 3 K o 



9^ 



9 y. 9 y 



(9) 



The solutions \|/ Q and ^J/ x of (9) should "be such 
that \|/ = y at infinity, according to (2) , that at f inity 
the straight lines y = ±l/2 H are streamlines ( H = tunnel 
width) , and lastly that streamline \|/ = 0 is a circle or 
nearly so. As solution ty 0 the one afforded by the con- 
ventional reflection method is taken: namely, 



H 2 TT 

7 - -ff- 



sin 



2-ny 
H 



Cosh^— cos^F 
a H 



(10) 



Formula (10) satisfies equation. M> 0 = 0 and com- 
plies exactly with the first two of the just quoted "bound- 
ary conditions. The contour ^ 0 = 0 is, Of course, not an 
exact cylinder, but the departure is slight as seen from 
the appended table, wherein r 0 , r 45 , and r 90 designate 

the values of radius vector r with respect to the points 
of the profile border incidental to * = 0° , 45°, and 90°. 



H/2 E 


r 0 /R 


r 45 /r 0 


r eo/ r o 


3 


0.9592 


0.99868 


0.99725 


5 


.9842 


.99979 


.99960 


10 


.9959 


.99999 


.99997 


20 


.9990 


1.00000 


1.00000 



The last two columns of the table reveal that the 
departure of the profile from circular form is practical- 
ly insignificant. (The departure of r 45 /r 0 , and r 90 /r 0 
from 1 is proportional to (B/H)?) Even in the extreme case 
of the profile spanning one-third of the tunnel width 
the thickness of the profile is only 0,27 percent less 
than the width. It is to be noted that E is not the 
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radius of the circle as the second column indicates. 

The solution of the second approximation reads 
^ = + a » and the prediction of \J/ X is predicated 
on the formulation of K 0 according to (9) and (10). 



K 0 = 1 - 



. „a a Cosh 2 -JI2 cos §JIZ. - 1 
4 S tt H H 



^Cosh ?IL X - cos 

1 

^Cosh - cos 



4 aV 4 



(11) 



Abbreviating these equations by means of 

H ' H H 

the last equation (9) assumes the form 



(12) 



9 u 2 3 v 2 



V 



= 2- sin v 



with 



L 



q s Cosh u q 4 - 2q 2 Sinh s u , q 6 Cosh u 

n „2 ' _ « —————— 



2 JJ B 2 N 3 

N = Cosh u - cos v 



8 N 



> (13) 



Then the solution of (13) is, expressed, in analogy 
with (10) in the form 



* x = - 



JL 

2tt 



gp (u) sin v 



(14) 



where H has the value given in (13); the g p (u) are 
functions of u only, and p is a suitably chosen whole 
number. The advantage of this formula is as follows: It 
at once satisfies the boundary conditions at the tunnel 
wall, v = ±tt for y = ±l/2 H according to (12); hence 
V x = 0 and i|r = \J/ Q = ±1/2 H. It is seen that formula (12) 
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is in fact practical and that >l> = ^ 0 = y for large x. 
Entering (14) in (13), H/2 it and, sin v afford abbrevi- 
ations, so that / 

p s j v 

V T (p " X) g P " V' + 8 P g*' slnh u £ 2p(p - 1) gpCosh u 1 

Li ~ H P ~" ~~ " r »p +1 "~~ " J 

q s Cosh u ^ q 4 - 2q 2 Sinh 2 u . q 6 Cosh u 

= — r~=5 + t + 3 — (15) 

2 B 2 2 N 8 H 

(the derivations with respect to u being indicated "by 
dashes) . 

Since v appears then solely in the function N = 
Cosh u - cos v, a comparison of the coefficients of equal 
powers of N on both sides of (15) gives the consecutive 
differential equations for the functions gp, and at the 
same time proves the possibility of formula (14). The 
denominator K' occurs only once in (I5)(that is, on the 
left side for p = 1); the correlated numerator must 
disappear: 

e/' = o 

As the consideration involves the flow past a profile 
symmetrical with both the X and Y axis, and hence all 
gp in x and also in u must be even functions. Therefore: 

gj = const = C 0 (16a) 

The correlated terms of N s give 

g 2 - g 2 » + 2 gl « Sinh u = jL q 2 Cosh u 

<> r g 2 n - g s = - \ q 2 Cosh u 

ge = Ci Cosh u - k q s xl sinh u (16b) 

Correspondingly the further functions g p follow at 
g 3 = C 2 Cosh 2 u + d + J- q 4 (16c) 
g 4 = C 3 Cosh 3 u + | D Cosh u (I6d) 

g 5 = C 4 Cosh 4 u + (2 C 3 + D) Cosh 2 u + — D fl6e) 

2 

g 6 = C 5 Cosh 5 u +^ | C 4 + | C 3 + | d) Cosh 3 u 

+ ( I °3 + 4~ ») Cosh u (I6f) 
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etc. With D abbreviating for 

D = C 8 + Ci + | q 4 * Cl6g) 

Quantities C ,-C 3 ,etc. are integration constants. 
For all values p |j 4. 

g p+1 » - p s g p+1 = 2 pgp'Sinh u - 2 p(p - l)g p Cosh u (17) 

As already stated "becomes zero for large value u; 

for large u the individual summands of (14) disappear as 
e" 11 , the still available integration constants are to he 
so defined that % - 0 disappears for the points of the 
original contour. The exact compliance with this require- 
ment affords an infinite series of functions gp; although 
it can he restricted to a finite series. For instance, 
with the selection of C 3 = - 3/2 D, C 4 = - l/2 C 3 - 1/4 D 
= l/2 D , and Cs = 0 , g 6 disappears identically, so that, 
according to (17) all suhsequent functions g 7 = g s = ••• 
can he put = 0. Formula (14) then terminates with a series 
of P = 5 terms with altogether three arhitrary constants 
C 0 , C x , C s , which are so determined that \|; x exactly dis- 
appears for the points of the profile of the first approx- 
imation (\j/ 0 = 0) related to * = 0° , 45°, and 90°. Then 
the profile of the second approximation agrees at least 
in eight points with that of the first. Exact agreement 
in more points can be attained hy increasing P. 

In the following P = 3 is deemed commensurate, the 
justification of this step being deferred till later. 
All the gp mu,st disappear identically for p l£ 4, which, 
according to (16) merely requires that D = 0 and further 
C3 = C4 = *» 0. Substitution of two new constants 

C 0 = c 0 q S ; Ci = c 1 q* (18a) 
for C 0 and affords 

0.-- (.1 ♦&)■«•- fc*' <!•»> 

from D = 0 according to (l6g) , with which values formula 
(14) gives, if the functions g- are taken from (16); 
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= - sin v 
2tt 



c q 2 c i<l* Cosh tt - i. q 2 u Sinn u 




(19) 



lext,' c. and c' are so defined that the square 
"brackets in (19) exactly disappear for the points of the 
profile edge (i$/ Q = 0) related to * = 0°, and * = 90°. 
For smaller values q series expansions, might be used to 
define u and v for these points according to (10) and 
then c 0 and c x according to (19).* 

The elementary calculation gives. 

1 



5 , 5 £ 
= — + — q £ 

12 48. 



+ — q 

144 



4 + q 6 

4320 



72 



+ -2— L_ q s 

1440 5670 



(20) 



In equations (19) and (20) is given. While c 0 

and c-l have been so defined that x becomes exactly 
zero for only four points of the original profile boundary, 
actually the thus obtained ty 1 disappears ident ically , as 
is readily proved, for all points of the original profile 
boundary^up to the terms vith q 4 inclusive. The profile 

+ i> 1 ct = 0 is therefore to be regarded with the same 
accuracy a.s the circle at the profile Q = 0, which proves 
the validity of the choice of P = S. 

Allowing the tunnel to become wider and wider, that 
is, let q approach 0, ^ x assumes the value 



sin 4 { 



_7_ + 1 lf_ 

12 2 J3 



1_ E 
12 r 4 



E sin 3 * 



V4 4 J/ 



(21) 



known from the theory of free flow of fluid past a circu- 
lar cylinder. 
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The velocity is defined according to (3). Since o 
relates, to £ - Uv equations (6) and (&) > restricted to 
quantities of the order of . a , afford 



P 



= l + i- U 0 - l) a 8 + 
hence (with the same degree of accuracy) 

W 0 



w x _ 
TT 



3y 



x L dx 2 . 9x 



3*o y - 1 . 

97 3 Sy~ 
9x J 



a' 



a 



s 



> 



(33) 



Disregarding, at this instance, the general expres- 
sion ^22) in favor of the narrowest - cross section, that 
is, the points x = 0 or u = 0, (9) and (19) give 

( w y^x=0 = 0 

+ C- 5. - al + _sl_ + i 6 1 a s 

V 4 13 120/ 4N s g 6N 3 J a ••• 



(23) 



with 



N 0 =l-cosv=3 sin 8 ' - 

2 



Two principal cases are considered. 

First, the maximum velocity at the tunnel wall (v v ) , 

where y = l/ 2 H , hence v - rr , tt 0 = 2; therefore * 

-#= 1 +■ | q S f [JJ q 2 + ^g* ... ] «■ . . ( 24a) 



second, the maximum speed at the profile edge (wp); for 
which (10) gives at * = 90° 



or 
w 



H 24 4 5760 



11 o 4 

q . . 



U 6 q + 180 q ' * * + L 6 



+ i§. q 3 + 4?_ 4 
72 * 432 q 



1 2 
• • • IOC • • • 



(24b) 
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In the 

n 4 



event that the tunnel is wide enough the terms 
with q 4 can be disregarded relative to those with q 2 , 
and H must then he put equal to the radius of the circle; 
likewise one may put a = M the Mach number M = U/a^ in 
undisturbed flow without introducing great error. 



Because (5) affords 



M = a 



a 



and the validity of the entire consideration holds only 
to about a = 0.4 since otherwise the maximum speed at 
the profile edge already reaches the value of the local 
velocity of. sound. Hence (24a) and (24b) may be expressed 
in the form . 



K 
U 



-&« 1 



4 



1 + ^" s 



s 



> 



(25) 



M = 0 furnishes the conditions in an incompressible 
fluid. The value 2 for the ratio w^/U is supplemented 
by a value depending upon the width of the tunnel. As a 
consequence of the compressibility, the value 2 and the 
supplementary value is raised in second approximation. 
The essential factor now is that the percentage increase 
in the additional amount is substantially greater than 

* ASK 2 

the rise of the basic amount 2, which is ™- M instead 

1*3 

7 a 

of ~r M , or more than six times as high. 

13 - • 

The "correction of the flow velocity" AU at the tun- 
nel wall is given by (26a): 



& 



(26a) 



The value of i D at the profile is, according to (25): 



w p = 2('U + 



Atr) ; [ 



1 + (M+AM) 



w^- ' .to be taken from (25). 



NACA Technical Memorandum No. 1018 



11 



Then 

(AU). 
U 

and consequently: 



(AU) p _ i 
(AU) K ~ 3 



(27) 



This means that no compressibility effect is perceptible 
in the ratio of the two AtJ in the employed approximation. 



Translation by J. Vanier, 
National Advisory Committee 
for Aeronautics. 
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